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Inelastic neutron scattering experiment has been performed on a new heavy fermion com- 
pound YbAgGe for a polycrystalline sample at zero magnetic field. A quasielastic scattering 
and a crystal field excitation at 12 meV was observed. Both of them are broadened with an 
intrinsic width of 0.9 meV at 1.5 K which increases to about 3 meV at 300 K. The temperature 
dependence follows To + AyT typically observed in Kondo systems. Crystal field parameters 
were also deduced. To explain the result that only one excitation is observed at 12 meV, it is 
concluded that the B22 term should be the main component. 
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A ternary rare-earth compound YbAgGe has been at- 
tracting interest as a new Yb-based heavy-fermion ma- 
terial. 1-5 One major aspect is the coexistence of Kondo 
effect and geometrical frustration in the fcagome-lattice- 
like triangular coordination of Yb atoms in ZiNiAl- 
type crystal structure (hexagonal space group P62m). 
Heavy electron mass is observed which reaches about 
400 mJ/mol K 2 at 2 K followed by successive magnetic 
phase transitions at quite low temperatures of 0.8 K and 
0.65 K. 1,2 The two step transition is considered to be a 
characteristic of the geometrical frustration as observed 
in isostructural YbPtln and YbRhSn. 6,7 Another aspect 
of importance is the magnetic field induced non-Fermi- 
liquid (NFL) behavior which is typically manifested in 
the — In T dependence of magnetic specific heat C mSLg /T 
in a wide temperature range from above 10 K to below 1 
K. 3 YbAgGe offers a good platform to study a crossover 
from magnetically ordered state with small moments to a 
Fermi-liquid (FL) state through a quantum critical point 
(QCP), as well as recently discovered NFL compound 
YbRh 2 Si 2 . 8 

Characteristics of Kondo effect and heavy fermion are 
typically manifested in the magnetic specific heat. 1 At 
low temperatures C mag /T increases with decreasing tem- 
perature in proportion to — In T from about 15 K to 3 K. 
Katoh et at analyzed the data with a single-ion Kondo 
model of J = 1/2 for 7k = 24 K. While the magnetic sus- 
ceptibility at high temperatures follows Curie- Weiss be- 
havior of well localized moment, x{T) for H \\ a exhibits 
a weak maximum around 4 K, which suggests disap- 
pearance of the local moment. Although this is reminis- 
cent of isostructural heavy fermion compound YbCuAl, 9 
YbAgGe is different in that a well defined crystal field 
(CF) anomaly is observed; C mag (T) exhibits a Schottky- 
type peak at around 60 K, which indicates the first 
excited CF level at around 110 K. 1 Our final goal in 
YbAgGe will be to understand the whole process from 
the local-moment formation at zero field to the FL state 
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at high fields through the QCP where the NFL state is 
observed. For this purpose, it is of essential importance 
to study correlation and fluctuation of the magnetic mo- 
ments by inelastic neutron scattering (INS). The aim of 
the present study, as a starting subject, is to investigate 
spin fluctuation associated with the Kondo effect and 
to determine the CF level scheme, using polycrystalline 
sample at zero magnetic field. 

INS experiment was performed using the LAM-D spec- 
trometer at the KENS pulsed neutron source. LAM-D 
is an inverted geometry spectrometer using crystal ana- 
lyzer mirrors. 10 The incident neutrons have a Maxwellian 
distribution of energy at room temperature and the scat- 
tered neutrons with an energy of 4.59 meV are selected 
by the (002) reflection of the pyrolytic graphite analyzer- 
crystals before being counted by the detector. The energy 
and momentum transfer are analyzed by the time of flight 
and the scattering angle. Two analyzer mirrors cover a 
range of low scattering angle around <f> = 35°; the other 
two cover high scattering angle around <fi = 85° which 
we use to estimate the contribution from phonon scat- 
tering and to check the scattering-vector dependence of 
magnetic scattering. 

A collection of small pieces of polycrystals with total 
mass 14.2 g was sealed in an aluminium container with 
helium exchange gas and attached to a liquid-helium 
cryostat. The details of the sample preparation and char- 
acterization are described in Ref. 1. The scattering inten- 
sity was put on an absolute scale in mb/meV/sr/Yb by 
comparing the intensity with that of the incoherent scat- 
tering from a vanadium standard sample. Contribution 
of phonon scattering was estimated from the data for 
<fi = 35° and (j> = 85° at high temperatures above 100 
K; at low temperatures it was estimated from the high 
temperature results by calculating the Bose factor. 

Figure 1 shows the magnetic part of the neutron scat- 
tering function S mag (<j), u) averaged for (j) = 35° and 85°. 
Two features are clearly observed; one is a quasielastic 
(QE) scattering around hw — meV and the other a 
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Fig. 1. Magnetic part of the neutron scattering function of 
YbAgGe per Yb ion. The data for 20 K, 50 K, and 100 K, are 
shifted by 15, 30, and 45, respectively. Solid lines are the fits 
with Lorentzian spectral functions. The dashed line is the mea- 
sured spectrum of a vanadium standard sample normalized to 
the intensity of the incoherent scattering of the YbAgGe sam- 
ple, representing the resolution function at meV. 

magnetic excitation around 12 meV. Both of them are 
much broader than the energy resolution. The 12 meV 
excitation can be ascribed to the crystal field (CF) exci- 
tation of the Yb 3+ ions. A scenario that the first excited 
doublet is located at 12 meV above the ground doublet is 
consistent with the Schottky-type specific heat anomaly 
around 60 K reported in Ref. 1 . No other significant peak 
structures were observed up to 60 meV. 

We have analyzed the magnetic excitation spectra with 

hw 2 
= l-e^(-WfeT)£ /iPi( ^ ;Ai) ' 

(1) 

where Pi represents a spectral function of the ith excita- 
tion with an energy transfer fku = A; and fi the strength 
of the excitation. Lorentzian spectral function was as- 
sumed here. Fitting results are shown by the solid lines 
in Fig. 1. With respect to the low energy excitation, the 
spectra above 20 K can be fitted well by the Lorentzian 
centered at hui — meV, which indicates that the excita- 
tion is quasielastic. Although there remains a possibility 
that the excitation is inelastic with A less than 1 meV, 
it was not possible to determine within our energy reso- 
lution. The excitation at 12 meV can also be fitted well 
by the Lorentzian at all temperatures. 

The temperature dependences of the spectral strengths 
fi roughly exhibit the behavior of the bulk magnetic 
susceptibility. The strength of the QE scattering de- 
creases with increasing temperature in accordance with 
the Curie term, while that of the 12 meV excitation do 
not vary with temperature very much since this corre- 



Fig. 2. Temperature dependence of the half width of the crystal 
field excitation and the quasielastic scattering. Solid and dashed 
lines represent fits with Tq + A*jT. 

sponds to the Van-Vleck term. The decrease in the ob- 
served intensity in Fig. 1 is due to the temperature effect 
in eq. (1). The sum of the two spectral strengths at 1.5 
K becomes ~ 160 mb, which corresponds to the static 
susceptibility of 0.042 emu/mol and is consistent with 
the bulk susceptibility reported in Ref. 1. 

Figure 2 shows the temperature dependence of the half 
width r of the Lorentzian spectral functions. The widths 
of both QE scattering and CF excitation become broad 
with increasing temperature. Although the CF excita- 
tion has broader width than the QE peak, this is due to 
the broader energy resolution at 12 meV than at meV; 
the two peaks have the same intrinsic widths. Even if we 
take into account the convolution effect by the resolu- 
tion function, the intrinsic widths of the QE peak will be 
almost the same as those in Fig. 2. 

Here, let us consider the CF level scheme. Since the 
eightfold multiplet of Yb 3+ (J — 7/2) splits into four 
Kramers doublets, three excitations are expected at the 
lowest temperature. However, only one excitation was ob- 
served experimentally. It is necessary to look for a level 
scheme in which only the excitation from the ground dou- 
blet to the first excited doublet has a significant tran- 
sition strength and the other excitations has negligibly 
small probability. 

Since the local symmetry of the Yb site is Civ when 
viewed from the a axis, all the CF parameters Bk q with 
even k and q up to six appear. Then, we first estimated 
the CF parameters from the electrostatic field of point 
charges in which we assumed negative charges for Ag 
and Ge, and a positive charge for Yb. As a result, we 
found that the B 22 term becomes the main source of 
the crystal field and the B 20 term modifies it slightly; 
higher order terms with k = 4 and 6 are negligibly small. 
This is a straightforward result of the crystal structure 
of YbAgGe. Looking from the c axis at the Yb ion at 
(a;, 0,0) on the 3/ site and the surrounding Ag and Ge 
atoms, we notice that the Yb ion is in a strong B22 type 
(= x 2 — y 2 type) crystal field. 11 The a axis becomes the 
principal axis of the CF anisotropy for this Yb ion. The 
anisotropy along the c axis, which is a natural result 
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Table I. Polycrystalline average of the squared matrix elements 
<7jjj}| 2 for the crystal field eigenstates of Yb 3 + described in 
the text. 
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Table II. Polycrystalline average of the transition strengths of the 



\i) -> \j) excitation at 1.5 K, 50 K, and 100 K. 
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(0,1) 
12.1 


(0,2) (0,3) (1,2) (1,3) 
19.7 27.8 7.6 15.8 
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T (K) 




transition strength (mb) 




1.5 K 
50 K 
100 K 


31.2 
27.3 
17.0 


0.51 0.055 0.00 0.00 
0.47 0.051 3.89 0.15 
0.33 0.038 8.6 0.40 


0.00 
0.44 
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of the hexagonal structure, induces the B 20 term. Since 
there is not an inversion symmetry, CF parameters of 
odd rank also exist; but these do not work as an actual 
crystal field, and the lowest order CF parameters are B 22 
and B 20 . 

Following the above observation, we have neglected the 
4th and 6th order terms and calculated the CF eigen- 
functions only with B 2 q and B 22 . First, with B 22 only, 
we found that only the transition strength between the 
ground and 1st excited doublets has a significantly large 
value. Next, we modified the cigenfunctions by slightly 
including B 20 so that the magnetization curve and the 
temperature dependence of the magnetic susceptibility 
can be reproduced as well as possible. Our model of the 
CF parameters are B 22 = -14.77 K and B 20 = 1.477 
K, which give 1st, 2nd, and 3rd excited levels at 140 K 
(= 12 meV), 228 K, and 323 K, respectively. 

The matrix elements of J±, the perpendicular com- 
ponent of J to the scattering vector, for thus obtained 
cigenfunctions are listed in Table I. It is noted that the 
direction of the principal axis of the CF anisotropy for 
the three Yb ions at (a;, 0,0), (0,2,0), and (— x,— x,0), are 
different from each other by 120°. In addition, the matrix 
elements were averaged over the direction of the scatter- 
ing vector to estimate the value for the polycrystalline 
sample. It can be seen from this table that the matrix el- 
ement for the (0, 1) transition, which represents the tran- 
sition between the ground state |0) and the 1st excited 
state |1), is much larger than those for the (0,2) and 
(0, 3) transitions. Then, at low temperatures where only 
the ground state is populated, it is expected that there 
appear only one sizable excitation for (0, 1). When the 
temperature is elevated so that the 1st excited level is 
populated, the excitation for the (1,2) transition is also 
expected. 

Calculated transition strengths at three different tem- 
peratures are listed in Table II in unit of mb. Corre- 
sponding energy positions are also listed. The calculated 
result at 1.5 K is consistent with the experimental re- 
sult that only one excitation is observed, and it can be 
concluded that the excitation corresponds to the (0,1) 



transition. At elevated temperatures above 100 K, an- 
other excitation corresponding to the (1,2) transition is 
expected to appear around 7.6 meV with its intensity 
about half as large as that of the (0, 1) transition. If this 
peak was observed experimentally, we could determine 
the energy level of the 2nd excited state; however, it was 
difficult to identify this excitation in the background of 
phonon scattering which dominates the observed inten- 
sity at high temperatures. More careful subtraction of 
phonon scattering using a reference sample is necessary 
to identify this peak. 

The experimental results and the analysis of INS sup- 
port the CF model in which the B 22 term is dominant. 
Figure 3 demonstrates the calculated magnetization and 
magnetic susceptibility for the CF parameters used in the 
analysis: B 22 = -14.77 K and B 20 = 1.477 K. Molecu- 
lar field is not considered; it just shifts l/x(T) along the 
vertical axis and changes the scale of magnetic field in 
the M(H) curves. Although the effective moment at high 
temperatures and the magnetic anisotropy along a and 
c axes are consistent with the experiment, 1/%(T) along 
the c axis at low temperatures deviates from the exper- 
imental results. Since the calculated ground state has 
strong magnetic anisotropy within the ab plane, x(T) 
for H || c levels off below temperatures where the 1st 
excited state is depopulated. However, the experimen- 
tal data follows the Curie- Weiss law down to the lowest 
temperature. Slight modification of B 20 cannot alter this 
behavior. Some mechanism which induce magnetic mo- 
ments along the c axis is necessary, which could be a 
hybridization with the conduction band or a ferromag- 
netic exchange in the Yb chain along the c axis. Con- 
cerning H || a, l/x{T) in calculation is not far from the 
experimental results if we include a molecular field. Cal- 
culated magnetization exhibits a saturation moment of 
about 2.5 /Ub, which also is roughly consistent with the 
experimental result that 2.0 /xb is induced at 15 T. 4 

These CF parameters correspond to the set No. 12 in 
Table I of Ref. 12 which can explain the bulk proper- 
ties and the INS results of an isostructural ErNiAl. 12 CF 
excitations in PrNiAl, NdNiAl, and ErCuAl, have also 
been measured, but the CF parameters have not been 
deduced. 13 It should be noted, in our analysis, that the 
sign of B 22 could not be determined since the calculated 
results little depend on the sign of B 22 in the present 
analysis. 

With respect to the spectral line shape, it is basi- 
cally described by a Lorentzian in the whole tempera- 
ture range of this study. The spectra in Fig. 1 may be 
regarded as single-site excitations. The QE line width 
follow roughly a r + A^/T law with intrinsic To = 0.9 
meV. The CF excitation is also broadened to the same 
width and exhibits the same temperature dependence. 
This behavior is commonly observed in Kondo systems 
in which the 4/ electrons are almost localized; the width 
at T = indicates a single-site spin fluctuation and is 
regarded as a measure of T K . 14 ' 15 The value of T = 0.9 
meV and T(300K) <~ 3 meV is slightly larger than the 
widths of YbBci 3 , YbAuCu 4 , YbPdCu 4 , and YbPtBi, 
which are in the Kondo regime with well localized mo- 
ments, 14, 16-18 but much smaller than those of YbCuAl 
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Fig. 3. Calculated magnetization curve at 5 K and temperature 
dependence of the inverse magnetic susceptibility for the crystal 
field model of YbAgGe described in the text. Experimental data 
in Rcf. 1 is shown by the circles (H \\ a) and squares (H | c). 



and YbAgCu 4 in the intermediate valence regime, 9 ' 14, 17 
Then, from the aspect of spin fluctuation, YbAgGe can 
be categorized in the Kondo regime and has a relatively 
high characteristic temperature 7k- 

The transition temperatures of 0.8 K and 0.65 K 
in YbAgGe are several times lower than the Neel 
temperatures of other isostructural heavy-fermion Yb- 
compounds with lower Tk: In =3.4 K and 1.4 K in 
YbPtln, 1.9 K and 1.7 K in YbRhSn, 3.5 K in YbPtSn, 
and 2.9 K in YbNiAl. 6,7, 19,20 In view of the fact that 
the Neel temperatures of RAgGe for R=Tb-Yb roughly 
scale with the de Gennes factors, 4 it is considered that 
the RKKY interaction gives Tn of a few Kelvin for the 
Yb compounds with ZiNiAl-type structure. Since a dis- 
tinctive magnetic structure of frustrated moments has 
been observed in YbNiAl which orders at 2.9 K, 21 it is 
supposed that the low transition temperature in YbAgGe 
is due to the Kondo effect rather than the geometrical 
frustration. 

At the lowest temperature of 1.5 K, the QE line shape 
deviates from the Lorentzian as can be seen in Fig. 1. 
This could be ascribed to emergence of spin correlation 
leading to a scattering-vector dependence of the QE in- 
tensity as observed, e.g., in CeCu6 and CeRu 2 Si2- 22,23 
This point will be an important subject in future exper- 
iments using a single crystalline sample. 

In conclusion, wc have measured inelastic neutron 
scattering spectra of a new heavy fcrmion compound 
YbAgGe for a polycrystalline sample at zero magnetic 
field. The spectra exhibit a quasielastic scattering and 
a crystal field excitation at 12 meV, both of which are 
broadened because of the spin fluctuation associated with 
Kondo effect. The temperature dependence of the line 
width follows a typical behavior of r + Ay/T observed in 
Kondo systems. The intrinsic width of 0.9 meV at 1.5 K 
and 3 meV at 300 K is relatively high in comparison with 



other Yb-based Kondo systems. Crystal field parameter 
was also deduced from the analysis of the observed ex- 
citations. Model calculation assuming the B22 term as 
the main component successfully explains the excitation 
spectra. This is related with the ZiNiAl-type structure 
in which the Yb site is in a strong field of anisotropy. 
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